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Abstract 

In this paper, a new simple method for accurate 
measurement of line frequencies is proposed. The method 
is based on storing the line frequencies to be measured in 
read-only memory ROM. For comparison, three 
alternative measurement methods are illustrated 

Keywords: Power line frequency, Read Only Memory, 
Power Control, Embedded Systems. 

1. Introduction 

There are different methods to measure line frequency of 
power systems. The most versatile of those are built 
around microprocessors. Microprocessors allow to add 
many functions beside measuring frequencies, such as 
adding an alarm, history and so on. Some of specialized 
frequency measurement appeared in the literature such as 
measurement under non stationary situations [6], and 
frequency mixing and sub-sampling [7]. Most of the time, 
however, it suffices to know the line frequency instantly 
and accurately. Our approach provides smart, cheap, fast 
and accurate method to measure the line frequency. As 
early ROM-based solutions allowed only small deviation 
from nominal frequency, recent advances in storage 
technology provides a huge accurate readings. Our 
approach provides instant measurement, because all what 
it takes is the propagation delay within the ROM. This, in 
particular, is important under instability conditions. We 
will present theoretical bases for the work, followed by a 
complete design and analysis of the circuit. At the end we 
will have conclusion and future research. 



2. Theory 

The operation of the proposed frequency meter is based 
on measuring the time period, T, of the line frequency to 
be found. The time period can be simply measured by 
counting the number of pulses A of a given clock 

frequency f in one time period. Therefore, the time 
period is given by 



It is obvious from this equation that the resolution in the 

1 

period measurement is — seconds. We will try in the 

f c 

following to device effective methods to carry the 
division operation of formula one in an effective and 
economic way. 

2.1. Processor-based Solution 

The processor-based frequency measurement is 
straightforward. Figure 3, presents a block diagram of a 
typical solution. The CPU counts number of pulses 
during one signal period. All operations including 
division and multiplication can be easily programmed, 
and the result in the desired format is produced by the 
CPU 10. Suitable and economic CPUs for this purpose 
do not have built-in division operations. Therefore, 
division is done through a sequence of shift and subtract 
operations, taking significant number of CPU cycles. 



The rest of the paper is organised as follows; Section 2 
presents the theoretical analysis of the subject which 
includes discussion about processor-based solution, the 
simplified frequency meter, ROM-based frequency meter 
basics and implementation. Section 3 formulates the 
conclusion. 



The CPU solution is versatile and it is easy to implement. 
The CPU can also be doing other useful work, such as 
issuing warning for certain frequency values, producing 
statistics of line frequency status, recording history of 
line frequency vibrations and others. 
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2.2. A simplified frequency meter 

In an earlier work, a simplified frequency meter has been 
devised .The division operation of Formula 1, has been 
approximated by a single addition operation. 



In Figure 4, The up/down counter with parallel load finds 
the difference, d , between the expected number clock 
pulses of the optimal frequency and the actual pulse 
count during one signal period. Hence, 



/ = 



i 



d 

T o ± V 

fc 



( 2 ) 



Using Taylor series expansion, and approximating by 
taking only the first two terms, the measured frequency 
becomes: 






T f 

1 0 Jc 



fo + a -d 



(3) 



It is not difficult to scale d by a certain constant. Finally, 
finding/ is turned into and add or subtract operation 
which can be done by a simple add/subtract circuit. 

2.3. ROM-based frequency meter 

For the measurement of a line frequency with known 
nominal value of f 0 Hz and frequency deviation of 

± D Hz, it is only needed to store a mapping table for 
the frequency range which extends from f 0 — D to 

f 0 + D Hz. The number of line frequencies M needed 

to be stored in the ROM is determined by the range of 
frequencies to be measured 2D and the resolution in the 
period of measurement: 

2 Df 

M = 7 7 r7 + 1 (4) 

JO U 



The value of the line frequencies /■ to be stored in the 



M locations of the ROM are given by 

, f e (fo+D) 

Ji f c + i(f 0 + D) 

With i = 0, 1,2,3, . .., M-l 



(5) 



The maximum line frequency f 0 + D Hz is stored in 



the first memory location having address 0 and the 
minimum line frequency f 0 — D Hz is stored in the last 



t memory location having the address M-l. 

By choosing the ROM size M to the power of 2 



( M =2”), the number of address lines will be n. The 
wordlength of each memory location is determined by the 
number of decimal digits representing the line frequency 
to be stored. For the simplicity of design, BCD format 
will be used to represent decimal digits. To store line 
frequency of m decimal digits, memory locations of 4m 
bits wordlength will be needed. 



Let us summarise the above discussion by designing a 
frequency meter for the line frequency of f 0 = 60 Hz 

with D = ±5 Hz deviation. Let us assume that 1024 line 
frequencies in the range from 55 to 65 Hz are to be stored 
in ROM. This requires a ROM with 1024 locations and 
wordlength of 4m bits. The number of address lines of 

ROM is 10, that is, (log 2 1024) From equation 2, the 

clock frequency is found to be 365722 Hz. This clock 
frequency gives a resolution in the time period 

1 

measurement of seconds. It is obvious from 

365722 

equation 3 that the resolution in the frequency 
measurement is not fixed and it has a worst-case 
resolution of 0.012 Hz. At 60 Hz. Therefore, in this case, 
four decimal digits are adequate to represent the line 
frequency. 

2.4. Implementation of ROM-based meter 

The purpose of line frequency meter can be easily 
implemented as shown in Fig 1. The main blocks in 
Figre 1 are 1 ) two one-shot multivibrators 2)?z-stage 
binary counter 3)ROM with ^-address lines and 4m - bit 
wordlength 4) m- 4-bit binary latches and 5) m seven- 
segment drivers and displays. The input signal whose 
frequency to be measured is applied to the conditioning 
circuit to produce TTL compatible signal. The TTL 
compatible signal is applied to a one-shot multivibrator 
(OS1) to generate pulses with very short duration (less 
than the time period of the clock frequency) with the 
same frequency as the input signal. The signal will be 
used to latch the output of the ROM. The output of OS1 
is applied to a second shot (OS2) to generate pulses with 
1 

duration equals to . This signal is used to clear 

f 0 +D 

and disable the counter as long as pulse is high. A stable 
high-frequency crystal clock whose frequency is 

determined from equation 2 given f 0 , M and D , is 

applied to the count-enable input of the counter. The 
timing diagram for the operation of the proposed 
frequency meter is shown in Fig 2. The counter counts 
only when the output of the OS2 is low and it is cleared 
and disabled when it is high. The number of counts in 
this time period represents the address of the line 
frequency which is stored in one of the M locations of the 
ROM. The line frequencies are stored in the ROM BCD 
format. The latched frequency is applied to the BCD-to- 
seven segment drivers followed by seven segment 
displays. The line frequencies is to be stored in ROM are 
shown in Table 1. It is obvious from the table that four 
decimal digits are needed to represent the line 
frequencies with an acceptable accuracy. 

3. Conclusion 

An accurate, low cost frequency meter for the 
measurement of line frequency is proposed. A ROM is 
used to store a range of line frequencies around the 
nominal line frequency. The number of line frequencies 
is determined by the acceptable deviation of the line 
frequency from the nominal value and the required 
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resolution. The proposed meter can be easily modified to 
measure other frequencies with known nominal values by 
changing ROM contents and clock frequency. Other bad 
conditions in line frequency, such as resisting harmonics 
or low/high inductance, or distortion will be left for 
future research. 
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Table 1. Example ROM contents 



Address 


Line 

Frequency 


0 


56.000 


1 


64.988 


2 


64.977 


3 


64.965 


467 


60.018 


468 


60.009 


469 


59.999 


470 


59.989 


1020 


55.025 


1021 


55.017 


1022 


55.008 


1023 


55.000 
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Figure 1. ROM-based line frequency meter 
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Figure 2. Timing diagram for the frequency meter. 
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Figure 3. CPU-based line frequency meter. 
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Figure 4. Simplified Frequency meter block diagram 
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Abstract 

Capacitor placement problem is a non-linear and non- 
differentiable mixed integer optimization problem 
with a set of equality and inequality constraints. Most 
conventional optimization techniques are incapable to 
solve this hard combinatorial problem. The radial 
distribution systems are unbalanced because of 
single-phase, two-phase and three-phase loads. Thus, 
load flow solution for balance radial distribution 
networks will not be sufficient to solve the 
unbalanced case and, hence, special treatment is 
required for solving such networks. Optimal capacitor 
placement can result in system loss reduction, power 
factor correction, voltage profile improvement, and 
feeder capacity release. In this paper, Particle Swarm 
Optimization (PSO) based approach is used to 
achieve optimal capacitor placement for balanced and 
unbalanced radial distribution systems. The size of 
capacitors can be determined by using particle swarm 
optimization algorithm. In this fitness function is to 
maximize the net savings. In order to obtain the 
fitness values, load flow calculations has to be carried 
and capacitor placement is almost certain that better 
loss reduction can be obtained. The proposed solution 
method employs to search for optimal locations, 
types, and sizes of capacitors to be placed and 
optimal numbers of switched capacitor banks at 
different load levels. The proposed approach has been 
implemented and tested on various test systems with 
promising results. 

Keywords: Capacitor placement, load flows, radial 
distribution system, Discrete Particle swarm 
optimization (DPSO) 

1. Introduction 

Due to high concentration of inductive loads in 
distribution system, energy losses are more. Shunt 
capacitors are widely used in distribution system 
consuming negative VAR, which counteracts some of 
the lagging components of inductive VAR at the 
point of installation. Thus, it modifies the 
characteristics of inductive load. Shunt capacitor 



results in to a number of benefits like improvement of 
power factor, reduction of power loss, improvement 
of voltage profile, improvement of voltage stability 
and system-released capacity. To achieve these 
benefits to the utmost extent under various operating 
constraints, distribution engineers are required to 
determine the optimal locations, types and sizes of 
capacitors to be placed and control settings of 
switched capacitors at different load levels. 

In the past lot of work has been carried in the area of 
reactive power compensation for radial distribution 
networks [1-3]. These methods are based on different 
non-linear programming techniques and provide only 
a local optimum solution with a less computation 
burden. Recently, researchers have attempted to 
obtain optimum values of shunt capacitors for radial 
distribution networks using Simulated Annealing and 
Genetic Algorithms (GA).Chiang have used the 
method of simulated annealing to obtain the optimum 
values of shunt capacitors for radial distribution 
networks. SundhaRajan and Pahwa [4] have used 
Genetic Algorithm for obtaining optimum value of 
capacitors. They [4] treated capacitors as constant 
reactive power load and no method is suggested to 
reduce CPU time. GA is capable of determining a 
near global solution with lesser computational burden 
than the simulated annealing method [4]. Capacitor 
placement problem is a non-linear and -differentiable 
mixed integer optimization problem with a set of 
equality and inequality constraints. Most 
conventional optimization techniques are incapable to 
solve this hard combinatorial problem. But particle 
swarm optimization (PSO) algorithm proposed by 
Eberhart and Kennedy is very suitable [5, 6]. The 
radial distribution systems are unbalanced because of 
single-phase, two-phase and three-phase loads. Thus, 
load flow solution for balance radial distribution 
networks will not be sufficient to solve the 
unbalanced case and, hence, special treatment is 
required for solving such networks. D.Thukaram, 
H.M.Wijekoon Banda and Jovitha Jerome proposed a 
robust three phase power flow algorithm for radial 
distribution systems [8]. Optimal capacitor placement 
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can result in system loss reduction, power factor 
correction, voltage profile improvement, and feeder 
capacity release. Xin-me Yu, Xin-yin Xiong, Yao-wu 
Wu proposed a PSO-based approach to optimal 
capacitor placement with harmonic distortion 
consideration [8]. 

In this paper, a modem Al-based method, particle 
swarm optimization is used to solve the capacitor 
placement with all the realistic problem formulation 
considerations. The size of capacitors can be 
determined by using particle swarm optimization 
algorithm for balanced and unbalanced radial 
distribution system. In this fitness function is to 
maximize the net savings. In order to obtain the 
fitness values, load flow calculations has to be carried 
and capacitor placement is almost certain that better 
loss reduction can be obtained. 

The paper is organized as follows: Section (2) 
focuses on load flow method for balanced and 
unbalanced radial distribution system. Section (3 & 
4) describes on development and performance of 
discrete particle swarm optimization (DPSO). Section 
(5) gives the applying DPSO for the capacitor 
placement problem. Section (6) gives the results and 
analysis & section (7) listing of the references used in 
the paper. 

2. Load Flow Method 

In this section load flow technique for balanced and 
unbalanced radial distribution systems are discussed. 

2.1. Balanced Radial Distribution System 

In any radial distribution system, the electrical 
equivalent of a branch- 1, which is connected between 
node 1 and 2 having a resistance r (1) and inductive 
reactance x (1) is shown in Figure. 1. From Figured 
current flowing branch- 1 is given by 



feihu 






| V (3)b% 

2 



r fl) + i s a 

j 0(3] 

Figure 1. Single line diagram of branch- 1 

! (1) = | V (l)| Z5 (l) " | V (2)| Z5 (2) /( r (l) + J X ( )) 



I (l) “ ( P (2) J’Q(2)) /V (2)* Z5 (2) 

From Equations. (1) and (2) 



P (2) 78(2) 

\ V (2)\ Z -<5(2) 



10 M 1) 



( 2 ) 



Separating real and imaginary parts, the real part is 

| V (1)|| V (2)| C0S ( S (1) “ S (2))= | V (2)| + { P (2) r (l) + Q(2) X (1)} P) 



and the imaginary part is 

|^(l)||^(2)| sin (^(1) _ ^(2))= {^(2) X (1) “ G( 2) r (l) } 

=> |^(2)| 4 + 2|V( 2 )| 2 (/Xl)^2) + *(1)2(2) - °-5|v(i)| 2 ) + 

(Ilf + *(lf)(^2) 2 + 2( 2 f) = 0 



Equation (5) has a straightforward solution and does 
not depend on the phase angle, which simplifies the 
problem formulation. In a distribution system, the 
voltage angle is not so important because the 
variation of voltage angle from the substation to the 
tail end of the distribution feeder is only few degrees. 
Note that from the two solution of |V( 2 )| 2 only the one 
considering the sign of the square root of the solution 
of the quadratic equation gives a realistic value. The 
same is applicable when solving for |V( 2 )|. Therefore 

from Equation (5), the solution of |v^ 2 j|can be 

written as 

| y (2)| = {{( r (l)^(2) + X (l)2(2) - 0 - 5 | y (l)| ) 2 



- ( r (i) 2 + x (i) 2 )i p (2) 2 + Q(i) 2 ) r (6) 

_ ( r (l)' P (2) + X (l)Q(2) ~ 0- 5 Kl)| ) I' 2 

In general 

K + 1)| = {{( r o) F o + i) + x o)Q(i + 1) - °- 5 K)l ) 2 - 
( r 0) 2 + x (i))i P (M) + 2(;+i) 2 ) } <7) 

- ( r 0) P ti + l) + X 0)2(; + 1) -°- 5 Ki)| ) } 



where, node no., i=l,2,...,nd 
branch no., j=l,2...,nd-l 
nd= total no. of nodes 

The active and reactive power losses in branch ‘j ’ are 
given by 



Ploss[j]= 



^{pVo + Q 2 (» i) 




} 



( 8 ) 



Qloss[j]= 






0+i) 




} 



(9) 



The total active and reactive power losses of 
the system are : 

nd - 1 

TPL= X Ploss[j] 

7=1 



( 10 ) 
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nd-l 

TQL= Yj Qloss[j] (ii) 

7=1 

Where Ploss[j], Qloss[j] =Active and reactive 
power losses of branch ‘j’ 

TPL, TQL = Total active and reactive 
power losses of the system 

2.2. Un-Balanced Radial Distribution System 



In many cases, it is observed that the radial 
distribution systems are unbalanced because of 
single-phase, two-phase and three-phase loads. Thus, 
load flow solution for balance radial distribution 
networks will not be sufficient to solve the 
unbalanced case and, hence, special treatment is 
required for solving such networks. To solve the 
recursive algebraic equations the following 
information is necessary: 

• Status of the feeder line, overloading of the 
conductor and feeder line currents. 



• Whether the system can maintain adequate 
voltage level for the remote loads. 

• The line losses in each segment. 

Circuit model is shown in Figure 2 for a 
three-phase, four-wire grounded star system. Line 
charges are neglected at the distribution voltage level. 
For this four- wire system, Carson’s equations lead to 
the development of an impedance matrix of 4x4 
dimension. This matrix is used to calculate conductor 
voltage drop as shown below. By using Krichhoff s 
voltage law, 



y 

i 

v bg 

i 

V c s 

i 




y ng _ y ng 



ze , 



„ ba „ bb „ be ^ bn 

ze ij ze tj ze . ze tj 



na nb nc nn 

ze a z,e v ze tJ ze # 





" I a ~ 




IJ 




I b 




IJ 




I c 




IJ 


1 


^ s 

1 



( 12 ) 



In a grounded neutral system, the voltages at neutral 
and ground are same, 

V ng _ y ng = o (13) 



-yag 


i 

i 




1 

• a 

1 




aa ab ac an 

zey ze tj ze tj ze tj 


1 

a 

1 


ybg 


_ybg 


= 


y; 


= 


ba bb be bn 

z^ ze tj zey zey 




yeg 


1 

1 




i 

a 

i 




ca cb cc cn 

\_zey zey ze tj ze tj J 


_n_ 



(14) 



Where the values of the impedance elements are 



computed using the following expression based upon 
Carson’s equation. 


ze:: xze H 

z aa = ze” J J 

y y 

ZC ij 




(15) 


Real and reactive power losses for phases 
can be written as: 


a , b and c 


SL-j = PL a + j QL“ = V : a x (l 




1;)' <■<>> 


SLy = PLy + j QLy = Vj b x (i 


gf-V-xl 


tj)' < 17 > 


SLy = PLy + j QL)j = V' x (i 


s)‘-v;xi 


(l‘J (18) 



Nogs i 



Node j 



w 


=; y X ) 




■ — 1 ► 

I 


f=T J T- 


m 


, r\ 


S>- J J 


w 


■ r ^ v — : 


Ptd \ 

? — 'T-? — J — J J V J 



Figure.2 Three phase four wire line model 



In general, the receiving end power at any phase a, 
of line between nodes i and j may be expressed as: 



p ii + J Qu 



Z i p k + J Qk ) 



k = index of all nodes fed 
through the line between 
nodes i and j 



+ 



(19) 



Z (c + j QK m ) 

index of all lines 
connected between 
nodes m and n fed 
through the line 
between nodes 
iand j 



by substituting Equation (13) in (12) 



Similarly for phases ‘b’ and ‘c’ respectively. 
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3. Development of Discrete Particle of the inequality, velocities were adjusted according 

Swarm Optimization t0 difference, per dimension, from best locations 



Particle swarm optimization as developed by the 
authors comprises concept, and paradigms can be 
implemented in a few lines of computer code. It 
requires only primitive mathematical operators, and 
is computationally inexpensive in terms of both 
memory requirements and speed. Early testing has 
found the implementation to be effective with several 
kinds of problems. Particle swarm optimization has 
also been demonstrated to perform well on genetic 
algorithm test functions. 

3.1 Simulating Social Behavior 

One motive for developing the simulation was to 
model human social behavior, which is of course not 
identical to fish schooling or bird flocking. One 
important difference is its abstractness. Birds and fish 
adjust their physical movement to avoid predators, 
seek food and mates, optimize environmental 
parameters such as temperature, etc. Humans adjust 
not only physical movement but cognitive or 
experiential variables as well. We do not usually walk 
in step and turn in unison (though some fascinating 
research in human conformity shows that we are 
capable of it); rather, we tend to adjust our beliefs and 
attitudes to conform to those of our social peers. 

3.2 Precursors: The Etiology of Particle Swarm 
Optimization 

The particle swarm optimization is probably best 
presented by explaining its conceptual development. 
As mentioned above, the algorithm began as a 
simulation of a simplified social milieu. Agents were 
thought of as collision-proof birds, and the original 
intent was to graphically simulate the graceful but 
unpredictable choreography of a bird flock. 

3.3 Nearest Neighbor Velocity Matching and 
Craziness 

A satisfying simulation was rather quickly written, 
which relied on two props nearest-neighbor velocity 
matching and “craziness.” A population of birds was 
randomly initialized with a position for each on a 
torus pixel grid and with X and Y velocities. At each 
iteration a loop in the program determined, for each 
agent (a more appropriate term than bird), which 
other agent was its nearest neighbor, then assigned 
that agent’s X and Y velocities to the agent in focus. 
Essentially this simple de created a synchrony of 
movement. 

3.4 Acceleration by Distance 

Though the algorithm worked well, there was 
something aesthetically displeasing and hard to 
understand about it. Velocity adjustments were based 
on a crude inequality test: If presentx > bestx, make it 
smaller; if presentx < bestx, make it bigger. Some 
experimentation revealed that further revising the 
algorithm made it easier to understand and improved 
its performance. Rather than simply testing the sign 



vx[ ][ ] = vx[ ][ ] + rand ( ) (2Q) 

xp incrementx (Pbestx [ ][ ] - presentx[ ] 

[]) 

(Note the parameters vx and presentx have two sets 
of brackets because they are new matrices of agents 
by dimensions; increment and bestx could also have a 
g instead of p at their beginnings.) 

4. Performance of Particle Swarm 
Optimization Using Inertia weights 

The following describes the position and velocity 
update equations with weight factors included. 



V id = W x V id + Ci x rand ( ) x ( P id - X id ) (?u 

+C 2 x RandQ x (P gd -X id ) 1 J 



Xid Xid + Vid 



( 22 ) 



Where C 1? C 2 are positive constants and called 
cognitive and social parameters 

Equation.(21) calculates a new velocity for each 
particle (potential solution) based on its previous 
velocity, the particle's location at which the best 
fitness so far has been achieved, and the population 
global (or local neighborhood, in the neighborhood 
version of the algorithm) location at which the best 
fitness so far has been achieved. Equation. (22) 
updates each particle's position in solution 
hyperspace. The two random numbers are 
independently generated. The use of the inertia 
weights, which typically decreases linearly from 
about 0.9 to 0.4 during a run, has provided improved 
performance in a number of applications. 

4.1 Parameter Selection in Particle Swarm 
Optimization 

Unlike many other computational intelligence 
techniques, the particle swarm optimizer has few 
parameters to tune. Many attempts have been made to 
improve the performance of originally developed 
PSO. Many parameters have been added to the 
originally developed PSO to modify or to improve 
the performance of the technique. A quick statistical 
experiment is used to fine tune these parameters for 
the class of constrained optimization problem 
considered. 

4.2 Swarm Size,P 

The swarm size or the number of individuals inside 
the population is determined by the integer parameter 
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‘P\ For very small values of P the possibility of 
being trapped in local optima is very likely. Larger 
population will increase the computation time 
requirements. 

4.3 Number of Iterations 

It is also a user defined parameter that determines the 
number of iterations for which algorithm has to run. 
Based upon computational experience a few hundred 
of iterations are usually sufficient to observe 
significant improvement in the solution providing 
that the initial solutions are feasible. 

4.4 Velocity of Particle 

One of the important factors that affect performance 
of the PSO, Specifically the speed of the convergence 
is the particle’s velocity of ‘flying’ inside the 
problem space. This parameter limits the steps taken 
by particles at every iteration. A small value can 
cause the particle to get trapped in local optima; on 
the other hand, a too large value can cause oscillation 
around a certain position. This problem of proper 
selection of velocity can be eliminated by using 
adaptive velocity to PSO’s position updating. The 
following equation represents the velocity update 
equation. 

V id = W x Vid + Q x rand() x (P ld -X id ) (23) 
+C 2 x Rand ( ) x ( P gd -X id ) 



C j( U j Hins, +k cj -^ (25) 

U s 

Where k inst is a constant representing a fixed 
installment cost of capacitor banks; k CJ is the purchase 
cost of one capacitor bank: k c j = k cf for fixed 
capacitor banks and k c j = k cs for switched capacitor 

banks with k cf < k cs . The cost function C ■ (u°j ) is a 

non-differentiable and piece-wise function as 
illustrated in Figure.4. Consequently, the capacitor 
placement problem is a non-linear mixed integer 
optimization problem with non-differentiable 
objective function that cannot be efficiently solved by 
conventional optimization techniques. 

▲ 



u 



o 



o 

3 

o 



No. of capacitors 



Where cl and c2 are positive constants, called 
cognitive and social parameters respectively, both are 
equal to 2 in general cases. ‘W’ is inertia weight 
factor; a large weight factor facilitates a global search 
while a small inertia weight facilitates a local search. 

5. Problem formulation for Capacitor 
Placement 

For a realistic distribution system with presence of 
nonlinear loads, the objective of optimal capacitor 
placement problem is to reduce the power loss and 
energy loss, and to minimize the cost of capacitor 
banks, while maintaining node rms voltages within 
prescribed values. The objective function can be 
expressed as 

minf(u°,u 1 )=k e £T i P i +£c j (u®) (24) 

i=l j=l 

The right-hand side of equation. (24) consists of two 
terms. The first term represents the energy loss cost 
and the second term represents the total capacitor cost 
which comprises the purchase cost and the 
installment cost. Practically, capacitors are grouped 
in banks of standard discrete capacities (e.g., a bank 
of 300 kVAr); therefore, capacitor sizes are regarded 
as integer multiples of the standard size of one bank. 
To describe realistic cost of capacitor placement, the 
cost function of capacitors at node j is 



Figure.4 The cost functions of capacitor banks. 



5.1. Illustrations of Candidate Node Identification 

The proposed node identification method is explained 
with 15- node system. After running the load flows 
for base case system, the active power loss is given 
by 60.34821 kW. By compensating reactive power 
injection at every node equal to local reactive load, 
one at a time and perform the load flows and 
calculate the active power loss in each case and also 
loss reductions. The results are tabulated in Table 1. 
The most suitable nodes for capacitor placement are 
chosen based on the condition PLI greater than a PLI 
tolerance value between ‘0’ and ‘1’. The tolerance 
value is selected by experimenting with different 
values in descending order starting from ‘ 1 ’ . The best 
value of the tolerance is the value which gives the 
highest profit and satisfying the system constraints. 
The Power Loss Indices (PLI) are calculated as 



_ (lossreductioifi] -minimunreduction) 
(maximum reduction- minimunreduction ] 



(26) 



5.1.1 Algorithm for Node Identification 

Following algorithm explains the methodology to 
identify the candidate nodes, which are more suitable 
for capacitor placement. 



Read radial distribution system data. 
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• Run the load flows and calculate the base 
case active power loss. 

• By compensating the reactive power 
injections (Q c ) at each node and run the load 
flows, to calculate the active power losses in 
each case. 

• Calculate the power loss reduction and 
power loss index using equation(26). 

• Select the candidate node whose 

PLI>Tolerance. 

• Stop. 



Table 1 Power loss reductions 



Node 

no 


Power loss after compensation Q c at 
each node (kW) 


Loss 

Reduction 

(kW) 


2 


58.948 


1.4005 


3 


56.982 


3.3665 


4 


53.078 


7.2704 


5 


57.876 


2.4719 


6 


57.376 


2.9721 


7 


58.41 


1.9379 


8 


54.54 


5.8081 


9 


55.311 


5.0368 


10 


57.587 


2.761 


11 


52.916 


7.4327 


12 


56.124 


4.2246 


13 


57.563 


2.7856 


14 


56.425 


3.9236 


15 


52.835 


7.5129 



Table 2 The power loss indices 15 node system 



Node 

no 


PLI 


2 


0.0000 


3 


0.3216 


4 


0.9603 


5 


0.1753 


6 


0.2571 


7 


0.0879 


8 


0.7211 


9 


0.5949 


10 


0.2226 


11 


0.9869 


12 


0.4620 


13 


0.2266 


14 


0.4128 


15 


1.0000 



6. Results and Analysis 

The proposed method is illustrated with two different 
cases as follows. 

Case-I: Illustrates the capacitor placement of 

balanced 15 -node radial distribution system. 

Case-II: Illustrates the capacitor placement of un- 
balanced 19-node radial distribution system. 

Case-I: The proposed algorithm is tested on 15-node 
radial distribution system. The line and load data are 
given in ref [9]. The power loss index (PLI) tolerance 
of 0.5 is chosen to get maximum profit as shown in 
Table. 3. The results of 15-node radial distribution 
system at different load levels before and after 
compensation by using fixed and switched capacitors 
are given in Table 4. The summary of test results is 
given in Table 5. The voltage profile before and after 
compensation are shown in Figure 4. From the 
results, it is observed that active power losses reduce 
from 60.34821 kW to 29.77601kW i.e., 50.6596 % 
loss reduction and minimum voltage improved from 
0.942389 p.u to 0.967561 p.u. thus voltage regulation 
improved from 5.761 1% to 3.2439%, i.e., 43.6930 % 
improved. 



Table. 3 Total profit in a year for different PLI tolerance values 



PLI 

tolerance 


No. of 
candidate 
nodes 


Candidate 

nodes 


Total 

profit 

(Rs) 


0.9 


3 


4, 11, 15 


900751.81 


0.7 


4 


4, 8, 11, 15 


964803.12 


0.5 


5 


4, 8, 9,11, 15 


979106.14 


0.3 


8 


3,4, 8, 9, 11, 
12, 14, 15 


947228.68 


0.1 


12 


3, 4, 5, 6, 8, 9, 
10, 11, 12, 13, 
14, 15 


343391.37 



Voltage profile for 1 5 Bus System 




Figure.4 Voltage profile for 15-node system before and after 
compensation for full load 
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Table.4 Test Results of 15-Node System at different load levels before and after Compensation 



Load 

level 


Capacitors (kVAR) 


Total Losses (kW) 


Minimum Voltage (p.u) 


Fixed 


Switched 


Before 

Compensation 


After 

Compensation 


Before 

Compensation 


After 

Compensation 


0.5 


Qci4)= o 


Qci 4 )= 300 


14.346233 


9.357631 


0.971951 


0.979091 


Qa8) = 0 


QasfO 


QavnO 




QcanrO 


QcmV=0 


Qcn5i = 0 


Qcasi = 0 


1.0 


Qci4)= 0 


Qci4)= 300 


60.348211 


29.776017 


0.942389 


0.967561 


Qcf8) = 0 


Qcm = 300 


QavnO 


Qcm = 300 


Qcan=0 


Qmn=300 


Qai5,=0 


Qcy.5,=0 


1.2 


Qavr 0 


Qcf41= 9 


88.774877 


44.877899 


0.930082 


0.955847 


Qci8) = 0 


Qa8)=300 


Qav> = 0 


Qcm = 300 


Qcmf= 0 


Qcair 300 


Qai5,=0 


Qcyi5i=300 



Table. 5 Summary of results of 15 node system before and after compensation for load level of ‘ V p.u 



Description 


15-node system 


33-node system 


Before 

Compensation 


After 

Compensation 


Before 

Compensation 


After 

Compensation 


Total reactive Power load (kVAR) 


1251.17 


51.17 


2300 


1100 


Reactive power loss 


59.3422 


28.9419 


135.239381 


96.375501 


Released reactive power Demand (kVAR) 




1230.4003 




1238.863880 


Real power losses (kW) 


60.348211 


29.776017 


202.706961 


143.725587 


Loss reduction (%) 




50.6596 




29.096866 


Real power Demand(kW) 


1286.7482 


1256.17601 


3917.706961 


3858.725587 


Released demand (kW) 




30.572194 




58.981374 


Feeder Capacity (kVA) 


1836.6173 


1258.7279 


4612.896993 


4039.935333 


Released feeder capacity (kVA) 




577.8893 




572.961661 


Min.voltage (p.u) 


0.942389 


0.967561 


0.913041 


0.925082 


Voltage regulation (%) 


5.7611 


3.2439 


8.6959 


7.4918 


Improvement of Voltage regulation (%) 




43.6930 




13.8468 


Power cost (Rs) 


1730786.69 


1033976.17 


5813635.63 


4302049.82 


Net savings (Rs) 




696810.52 




1511585.81 



Case-II: The proposed algorithm for unbalanced 
radial distribution system is tested on 19-node 
radial distribution system. The line and load data of 
19 node radial distribution system are given_ref [7]. 
The summary of test results of 19-node radial 
distribution system is given in Table. 6. From the 
results, it is observed that total active 
power losses reduce from 13.473 kW to 1 1.257 kW 
i.e., 16.443 % loss reduction and minimum voltage 
of phase A, phase B and phase C 
improved from 0.95159, 0.94975, 0.95046 p.u to 
0.95751, 0.95572, 0.95641 p.u. respectively. The 
voltage profile for phase A, phase B and phase C 
before and after compensation are shown in 
figure. 5., figure. 6 and figure. 7. respectively. 

7. Conclusion 

A Discrete Particle Swarm Optimization for 
Optimal Capacitor Placement in Radial 
Distribution System has been proposed in this 



paper. With hill considerations of different load 
levels, and practical aspects of fixed or switched 
capacitor banks, the target problem is reformulated 
by a comprehensive objective function and a set of 
equality and inequality constraints. The proposed 
solution method employs PSO to search for optimal 
locations, types, and sizes of capacitors to be 
placed and optimal numbers of switched capacitor 
banks at different load levels. 

An effective approach for optimum location of 
capacitor in radial distribution system has been 
proposed. From the results, several important 
observations can be concluded as follows. 

• The power losses of distribution system can be 
efficiently reduced by proper placement of 
capacitor. 

• In addition of power loss reduction, the voltage 
profile can be improved as well by the 
proposed method. 
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Table.6 Summary of test result of 19-node system before and 
after compensation 



Description 


1 9-node system 


Before 

Compen 

sation 


After 

Compen 

sation 


Qc required (kVAR) 


90.000 




Total reactive 
Power load (kVAR) 


183.066 


101.113 


Released reactive Power 
(kVAR) 




81.953 


Min. Voltage (p.u) 


PhaseA 


0.95159 


0.95751 


PhaseB 


0.94976 


0.95572 


PhaseC 


0.95046 


0.9565 


Voltage 
Regulation (%) 


PhaseA 


4.8408 


4.2486 


PhaseB 


5.0242 


4.44281 


PhaseC 


4.9534 


4.3591 


Improvement of 
Voltage regulation 


Phase A 




12.2335 


PhaseB 




11.5717 


PhaseC 




11.9978 


Total losses (kW) 


13.473 


11.258 


Loss reduction (%) 




16.443 


Demand (kW) 


379.413 


377.198 


Released demand (kW) 




2.215 


Feeder capacity (kVA) 


421.268 


390.515 


Released Feeder Capacity(kVA) 




30.753 


Net savings (Rs) 




38012.06 



Voltage profile for Phase A 




0.955 

0.95 1 ' 1 ' 1 1 1 1 ' 1 ' ' 1 ' ' ' ' 1 ' 

0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 

Node number 

Figure. 5. Voltage profile for Phase A 



Voltage profile for Phase B 




Node number 



Figure. 6. Voltage profile for Phase B 



Voltage profile for Phase C 




Figure. 7. Voltage profile for Phase C 
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Abstract 

Discrete sliding mode control (DSCM) is one of the 
best techniques in analyzing the dynamics of 
nonlinear systems. This paper presents an improved 
discrete quasi- sliding mode control strategy in which 
the control input for the next step is determined by 
considering the dynamics of the current disturbance, 
which improves the robustness of the system 
dynamics with relatively smaller control effort. In 
fact, the existing techniques fail to deliver 
satisfactory response for parameter variations and 
time- varying external disturbances, where as the 
complete robustness is guaranteed in the proposed 
technique. This improved strategy is applied to two 
different systems under various environmental 
conditions and the results obtained demonstrate its 
effectiveness and improved efficiency. 

Key Words : Variable -structure system, discrete time 
system, quasi sliding-mode, state-space. 

1. Introduction 

Variable structure control (VSC) is a general 
approach for the design and control of a class of 
nonlinear systems. VSC is implemented along with 
the stochastic control techniques like fuzzy logic, 
neural networks to obtain better dynamics [1]. In 
variable structure control the system is allowed to 
vary its structure by properly and deliberately 
changing the sign and / or magnitude of the input, 
forcing discontinuities in the input with respect to 
time. The plane which separates these different 
structures is called as switching plane or switching 
surface. This discontinuous input makes the phase 
trajectory of the system to undergo two modes, 
reaching mode and sliding mode. In reaching mode, 
the system phase trajectory, starting from anywhere 
on the phase plane moves toward a switching plane 
and reaches it in finite time. This is followed by 
sliding mode in which the phase trajectory 



asymptotically tends to the origin of the phase plane 
[6], [8], [9]. The switching surface decides the closed 
loop dynamics of the system which are at the 
designer’s choice and hence, also known as sliding 
mode control (SMC). The main advantage of SMC is 
its insensitivity to parameter variations, external 
disturbances and modeling errors [6], [9]. 

With the advent of digital computers and its 
widespread use in control systems, considerable 
efforts have been put in the study of discrete time 
VSC/SMC techniques, called as discrete sliding mode 
control techniques (DSMC) [2], [3], [6]. In discrete 
sliding mode the control input is constant over 
sampling periods. Hence, when the states reach the 
switching surface, the subsequent control would be 
unable to keep the states to be confined to the surface. 
This leads the system to undergo only quasi- sliding 
mode, i.e., the system states would approach the 
sliding surface but would generally be unable to stay 
on it [6]. Thus, DSMC does not posses the invariance 
property found in continuous time sliding mode. 

Gao introduced ‘ reaching law ’ method to design the 
controller for continuous-time VSC [4] and extended 
the same for discrete time counterpart [3]. This 
approach is found satisfactory when compared to the 
other methods proposed in [5]. The switching 
function s(/c) is effectively controlled to meet the 
required dynamics and also to satisfy the constraints 
of DSMC. This is followed by the derivation of the 
control law in conjunction with the known plant 
model and parameter variations. 

The technique proposed by Gao is simpler and 
directly deals with the reaching process and makes it 
easy to obtain the control law [4]. However, 
chattering in the steady state is a major drawback, 
which is due to the discontinuous switching control 
applied to the plant, which excites the unmodelled 
high frequency dynamics of the system. 
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Bartoszewicz proposed a state-feedback-based 
control law for uncertain systems with bounded 
uncertainties that guarantee’s discrete sliding mode 
[2]. This is an improved version of Gao’s method and 
ensures finite time convergence of phase trajectories 
on the sliding plane without chattering. However, this 
method is effective only for systems with time 
invariant external disturbances. For systems with 
parameter uncertainties and time varying external 
disturbances this technique fails to deliver 
satisfactory results as the system becomes unstable or 
stable with steady-state errors. 

By redefining the priori known function defined in 
[2], it is found that this technique requires lesser 
control effort, but still unsuitable for systems 
influenced by parameter variations and external 
disturbances. In this paper a redefined discrete quasi- 
sliding mode control strategy based on 
Bartoszewicz’ s control law for unbounded and time 
varying uncertainties is proposed. The modified 
technique not only improves the robustness of the 
system dynamics, but also requires smaller control 
effort compared to the previous strategies, while 
retaining the finite time convergence of the state 
trajectory to the sliding surface. 

The paper has been organized into six Sections. 
Section 2 gives an overview of DSMC and reaching 
condition details. Bartoszewicz technique is briefly 
discussed in Section 3 along with the redefined 
version of it. Section 4 explains the new improved 
quasi- sliding-mode control strategy for an uncertain 
system with unbounded uncertainties and the details 
of the two numerical examples considered in this 
paper. The responses of the example systems under 
different conditions and the improvements in the 
system dynamics and efficiency are discussed in 
Section 5. Section 6 provides the conclusion and 
future work. 

2. Discrete Quasi-Sliding Mode and 
Reaching Conditions 

A. Preliminaries: 

Consider a discrete-time system represented by, 

x(k + 1) = Ax(k) + AAx(k ) + bu(k ) 

+ /(*) ( 1 ) 

y(/c) = h T x(k ) 

where x is the n X 1 state vector, A is an n x n 
system matrix, b and h are input and output vectors 
of appropriate dimensions, u is the system input, and 
y is the system output. The pair (A, b) must be 
controllable. Further, the n x n parameter variation 
matrix AA and the n x 1 external disturbance vector 
/ satisfy the matching conditions [6], 

AA = bA A - a row vector 



f = bf /- scalar 
Define the switching function 

s(k ) = c T x(k ) (2) 



with vector c such that c T b ^ 0 and the resulting 
quasi-sliding motion is stable. Disturbances and 
parameter variations are bounded, so that the 
following relation holds: 

d L < d(k) = c T AAx(k) + c T f(k ) < d u (3) 

where the lower bound d t and the upper bound d u are 
known constants and d(k ) is the disturbance vector. 
Further more define d 0 = 0.5 (d t + d u ) and 8 d = 
0.5 (d u — di), where d 0 the average is value of 

d(k) and 8 d is the maximum possible deviation as 
shown in Figure. 1. The quasi-sliding mode is defined 
as the motion such that \s(k)\ < e , where the positive 
constant e is called the quasi- sliding-mode 

bandwidth. 



d\ d 0 d u 




Figure. 1. Deviation of d(k) 

B. Reaching condition and reaching mode: 

The condition under which the system states starting 
from any initial state, move towards the sliding 
surface and reach it in finite time, called as reaching 
condition or reaching law. The system trajectory 
under the reaching condition is called the reaching 
mode or reaching phase. So, under ideal conditions 
there exists a finite time t e such that for all t > t e , 
the sliding function s(x) = 0. 

In continuous time the reaching law is a differential 
figure which specifies the dynamics of a switching 
function s(%). A simple reaching law may be written 
as [4] 

ds ( ( w (4) 

— = - q(sgn(s )) 

In the above condition the parameter q decides the 
rate of convergence of s(x) towards the sliding 
surface and the magnitude of chattering in the sliding 
mode. Hence, by proper choice of the parameter q in 
(4), the dynamic quality of VSC system in the 
reaching mode can be controlled. 

In case of discrete mode the condition given in (4) is 
written as 



